Quality control by a mobile molecular workshop: quality versus quantity 
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Ribosome is a molecular machine that moves on a mRNA track while, simultaneously, polymer- 
izing a protein using the mRNA also as the corresponding template. We define, and analytically 
calculate, two different measures of the efficiency of this machine. However, we arugue that its 
performance is evaluated better in terms of the translational fidelity and the speed with which it 
polymerizes a protein. We define both these quantities and calculate these analytically. Fidelity is 
a measure of the quality of the products while the total quantity of products synthesized in a given 
interval depends on the speed of polymerization. We show that for synthesizing a large quantity of 
proteins, it is not necessary to sacrifice the quality. We also explore the effects of the quality control 
mechanism on the strength of mechano-chemical coupling. We suggest experiments for testing some 
of the ideas presented here. 

PACS numbers: 87.16.Ad, 87.16.Nn, 87.10.Mn 



I. INTRODUCTION 

For cyclic machines with finite cycle time, the effi- 
ciency at maximum power output [l|-|3[ is a suitable indi- 
cator of its performance For molecular motors, the 
Stokes efficiency @ is an alternative measure of per- 
formance. However, not all machines are designed for the 
sole purpose of performing mechanical work. In this pa- 
per we consider a class of machines whose main function 
is to synthesize a hetero-polymer, subunit by subunit, 
using another hetero-polymer as the corresponding tem- 
plate. 

The conceptual framework that we develop here is gen- 
erally applicable to all machines which drive template- 
dictated polymerization. But, for the sake of concrete- 
ness, we formulate the theory here for a specific machine, 
namely, the ribosome ff\; it polymerizes a protein using 
a messenger RNA (mRNA) as the corresponding tem- 
plate and the process is referred to as translation (of ge- 
netic message). The subunits of a mRNA are nucleotides, 
whereas those of a protein are amino acids. The cor- 
rect sequence of amino acids to be selected by the ribo- 
some is dictated by the corresponding sequence of the 
codons (triplets of nucleotides). The ribosome also uses 
the mRNA template as the track for its movement. In 
each step, the ribosome moves forward by one codon on 
its track, and the protein gets elongated by one amino 
acid. Ribosome has a quality control system to minimize 
translational error by rejecting incorrect incoming amino 
acid subunits. 

For a ribosome, the average speed V of polymeriza- 
tion of a protein and the fidelity <f> of translation, rather 
than efficiency and power output, are the primary indi- 
cators of its performance. We define these quantities, as 
well as a few others, quantitatively and calculate these 
analytically. Contrary to naive expectations, we show 
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FIG. 1: Pictorial depiction of the full chemo-mechanical cycle 
of a single ribosome (see the text for details). 



that a higher rate of elongation of the protein does not 
necessarily sacrifice the translational fidelity. 



II. THE MODEL 

A ribosome consists of two interconnected parts called 
the large and the small subunits. The small subunit binds 
with the mRNA track and decodes the genetic message 
of the codon, whereas the polymerization of the pro- 
tein takes place in the large subunit. The operations 
of the two subunits are coordinated by a class of adapter 
molecules, called transfer RNA (tRNA). One end of a 
tRNA helps in the decoding process by matching its anti- 
codon with the codon on the mRNA, while its other end 
carries an amino acid subunit; in this form the complex 
is called an amino-acyl tRNA (aa-tRNA). 

The three main steps of each cycle of a ribosome are 
as follows: (i) selection of the cognate aa-tRNA, (ii) 
formation of the peptide bond between the amino acid 
brought in by the selected aa-tRNA and the elongating 
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protein, and (iii) translocation of the ribosome by one 
codon. However, some of these steps consist of impor- 
tant sub-steps. Moreover, the branching of paths from 
a well defined mechano-chemical state can give rise to 
the possibility of more than one cyclic pathway for a ri- 
bosome in a given cycle. Neither of these features was 
captured by the earlier models [1, 1^ . The complete set 
of states and the pathways for a single ribosome in our 
model are shown in Fig[TJ This model may be regarded 
as an extension of the Fisher-Kolomeisky generic model 
for molecular motors [Toj . 

The selection of aa-tRNA actually consists of two sub- 
steps. In the first sub-step, among the tRNAs, which 
arrive at the rate uJa, non-cognate ones are rejected, 
at the rate w^i, because of codon- ant ico don mismatch. 
The second sub-step implements a kinetic proofreading 
mechanism for screening out the near-cognate tRNAs; 
this sub-step is irreversible and involves hydrolysis of 
a Guanosine triphosphate (GTP) molecule (at the rate 
cj/ii). Near-cognate tRNAs are rejected at the rate uJr2- 
But, if the tRNA is cognate, the protein gets elongated 
by the addition of the corresponding amino acid (at the 
rate oOp). However, occasionally a near-cognate tRNA 
escapes potential rejection and the ribosome erroneously 
incorporates it into the protein (at a rate flp). The 
subsequent translocation step actually consists of impor- 
tant sub-steps. The first sub-step is a reversible relative 
(Brownian) rotation of the two subunits with respect to 
each other (with rates Wf, / and uJbr along the correct path 
and with rates ilby and flbr along the wrong path). In 
the state labeled by 4 (and 4*), the tRNAs are in the 
so-called "hybrid" configuration 0] , the details of which 
are not required for our purpose here. The second sub- 
step, which is irreversible, is driven by the hydrolysis of 
a GTP molecule; this sub-step leads to the coordinated 
movement of the tRNAs within the ribosome and the 
forward stepping of the ribosome by one codon on its 
mRNA track (at the rates uJh2 and r2/,.2 along the correct 
and wrong paths, respectively). Further detailed identi- 
fication of the states labeled in FigUby the integers Hil 
is not needed for our purpose here. 

Thus, following the selection of a cognate aa-tRNA, 
the correct pathway isl— >2— >3^'4— >5— In 
contrast, if a non-cognate aa-tRNA is picked up, the most 
probable pathway is 1 2 — !• 1. However, if the aa- 
tRNA is near cognate, then the pathway could be cither 
1 ^ 2 ^ 3 ^> 1 (successful kinetic proofreading) or 
1— >2— >3^>4*— >5*— >1 (incorporation of a wrong 
amino acid). 

In principle, the model does not necessarily require any 
relation among the rate constants. However, throughout 
this paper we assume that 



We define the fraction 
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as a measure of translational fidelity. The error ratio 



^p/{ujp + flp) 
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is the fraction of wrong 



amino acids incorporated in a protein. We also define 
the rejection factor as 
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Note that the average velocity V oi a ribosome is also 
identical to the average rate of elongation of the protein 
that it polymerizes. A higher V results in a larger "quan- 
tity" of the protein after a given interval. On the other 
hand, the parameter (p characterizes the "quality" of the 
final product while TZ characterizes the "effort" of the 
quality control system in screening out the non-cognate 
tRNA (including near-cognate tRNA). 



III. RESULTS 

Suppose that 'P^{t) is the probability of finding the ri- 
bosome in the "chemical" state /i at time t, irrespective 
of its position on the mRNA track. The obvious normal- 
ization condition is X]^=i 'P^ + ^4 + ^5 = 1- Solving the 
master equations for V^{t) in the steady-state, we get the 
average velocity 



V = ic{uJh2P5 + ^h2Vl) = IcKJf + ^ j 



(4) 



where Ic is the length of a codon and 
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Note that in the special case VLp = 0, different terms 
iTgyy are the average time spent by the ribosome in 
different mechano-chemical states. 



of K 
iffer 

The direct transition 3^1 gives rise to "slippage" 
[l^ . Therefore, the coupling between the chemical input 
and the mechanical output is "loose" 113. The ratio of 
the mechanical flux to chemical flux [14| can be taken 
as a measure of the strength of the mechano-chemical 
coupling k: 



a constant, because the more stringent the proofreading 
is, the fewer will be the wrong proteins produced. 
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FIG. 2: The average velocity V of a ribosome (in the units 
of "codons per second"), i.e., average rate of elongation of a 
protein (in the units of "amino acids per second"), is plotted 
against (a) the rejection factor 7?,, (for six different fixed values 
of the parameter cjp), (b) the fidelity i\> (for four different fixed 
values of the parameter 01^2), and (c) the rejection factor TL 
(for four different fixed values of the fidelity 0) . 



SO that K = 1 in the hmit Wr2 — 0. Similar quantitative 
measures of mechano-chemical coupling have been intro- 
duced earlier [Tsj in the general context of motors which 
can execute "futile" cycles of hydrolysis of nucleotide tri- 
phosphates (including GTP). In the steady state of our 
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FIG. 3: The mechano-chemical coupling strength k is plotted 
against ijjr2 for four different fixed values of the parameter Wp. 




FIG. 4: (Color online) Thermodynamic efficiency r]T and 
Stokes efficiency r\s are plotted against the fidelity <j!>. 



model, 



2(wp -I- ilp) _ 2(a;p -f C - 01^2) 



2(a;p -f rip) + 'jJr2 2(a;p + C) - U}r2 
The input power consumed in proofreading is 

UJr2 



P? = 



Afj. 



(7) 



(8) 



whereas the input power consumed in protein synthesis 
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where A/i is the free energy released by the hydrolysis of 
a single GTP molecule. Thus, the total power input is 



p^^ ^{ -r2 + 2u;p + 2np ^ ^^^^^^ 



(10) 



We define the thermodynamic efficiency rjT and the 
Stokes efficiency rjs by the relations 



FV , -fV^ 
VT = and ris = — — 



(11) 
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where F is the load force and 7 V is the phenomeno- 
logical form of the viscous drag on the ribosome. One 
of the Umitations [3] of definitions (fTTj) is that it impUes 
vanishing thermodynamic efficiency as the external load 
force vanishes although, in reality, the motor continues 
to work against the viscous drag. 77s is the appropri- 
ate measure of the efficiency of a molecular motor when 
F vanishes. From definitions pT|) -qs = {jV/F)riT, and 
using expression (jj]) for we get 



{uJr2 + '^OJp + 2ilp)A/X 

Moreover, we assume the standard form 10] of in- 
dependence of uJh2{F) and $7^2 (-F): the values of these 
rate constants in the absence of F are multiplied by the 
factor exp{—S F£c/kBT) where < S < 1. For plotting 
the graphs, we have used the parameter values F — IpN, 
A/x = lO/csT, 4 = 0.9nm, S = 0.5, uja = 25.0s-\ 
Whi = 25.0s-\ Wft2(0) = 25.0s-\ Wf,,. = uj^f = lO.Os"!, 
ujp = 25.0s~^, and 7 = 60pN.s.m~^. The parameters 
LUri and ujr2 have been varied over a wide range always 
keeping ujr2 + = 20s~^. 

The average rate V of elongation of a protein is plotted 
against the rejection factor TZ and fidelity (f> in FigsIH^a) 
and (b), respectively (note the unit used for plotting 
V). The increase of ujr2 increases the rejection factor TZ. 
Therefore, naively, one would expect V to decrease with 
increasing TZ. However, because of constraint ([IJ, the 
increase in ujr2 is compensated by decreasing flp. There- 
fore, instead of decreasing, V can increase with TZ pro- 
vided ujp is sufficiently large (see Figl^fa)). 

Moreover, contrary to naive expectation, V increases 
with increasing (f> as long as ujr2 is kept fixed (see 
FiglSJb)). However, for a given (f>, V decreases with in- 
creasing ujr2', this decrease is caused by the increasing 
frequency of the futile cycles in which a molecule of GTP 
is hydrolyzed in kinetic proofreading and the aa-tRNA is 
rejected. In order to emphasize the interplay of TZ and 
(j) directly, in Figl5](c), we plot V against TZ for several 
different values of the parameter (f>. Clearly, for a given 
(j), V decreases monotonically with TZ. Similar monotonic 
decrease in V with increasing u;r2 demonstrates that the 
"slippage" caused by the kinetic proofreading weakens 
the mechano-cheniical coupling k (see Fig|31). 

The variation of the efficiencies rjT and r/5 with is 
shown in the three-dimensional plot of FigHl Taking 



cross sections of this diagram parallel to the rjT — ris plane 
for several different constant values of </) (not shown in 
any figure) we find that t]t increases monotonically with 
rjs, but the increase is sublinear. This is a consequence of 
the fact that although cj) depends only on the ratio Qp/ujp, 
V depends separately also on iVp through its dependence 
on Keff. 

IV. SUMMARY AND CONCLUSION 

Intracellular machines which carry out template- 
dictated polymerization of macromolecules can be re- 
garded as molecular motors [7, 16]. Obviously, as we have 
shown here, appropriate measures of efficiency can be de- 
fined and calculated for such motors. However, the per- 
formance of these machines is evaluated better in terms 
of fidelity and rate of polymerization, rather than effi- 
ciency and power output. Using our theoretical model 
for the kinetics of translation, we show that, contrary to 
the widespread belief, a higher rate of polymerization of 
a protein does not necessarily compromise the transla- 
tional fidelity. We have also analyzed the interplay of 
the "futile cycles" arising from kinetic proofreading and 
the consequent "looseness" of the mechano-chemical cou- 
pling. 

For carrying out the calculations analytically, we ig- 
nored the possible variations of the rate constants from 
one codon to another. Therefore, our prediction can be 
tested in-vitro using an artificially synthesized sequence- 
homogeneous mRNA (i.e., mRNA strand whose codons 
are all identical, except for the start and stop codons) . In 
the surrounding medium only two species of amino acid 
subunits should be available, one of these two species is 
cognate whereas the other is not. The rate constants ujri, 
<^r2, ^p, ^bf, ^br and fth2 cau be varied by replacing the 
non-cognate amino acid subunits with another distinct 
species which is also non-cognate. 

The performance of all the intracellular machines of 
template-dictated polymerization can be characterized 
and evaluated within the general conceptual framework 
developed here. However, the results for polymerases 
[l7| , which polymerize DNA and RNA, would differ from 
those of ribosomes because of the differences between 
their respective mechano-chemical cycles. 

We thank Joachim Frank for constructive criticism of an 
earlier draft of the manuscript. 
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